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INTRODUCTION 


For  the  past  two  years  the  collaborative  effort  between  the  University  of 
Poitiers  and  ARRADCOM  has  been  concentrating  on  the  physics  of  the  shock  front 
rise  time  (structure)  in  liquids.  Liquid  water  has  been  the  primary  vehicle  for 
that  study. 

Recent  success  in  obtaining  optical  reflectivity  data  for  liquid  nitrometh- 
ane  and  in  writing  a  more  complete  reflectivity  theory  (rereflection  has  now  been 
included)  have  shifted  the  program's  emphasis  from  the  shock  front  rise  time  to 
the  systematic  differences  observed  in  the  shock  front  reflectivity  of  water  and 
nitromethane  (for  assumed  comparable  rise  times).  This  is  the  first  report  of 
the  more  complete  theory  and  those  systematic  differences. 

Any  systematic  difference  between  the  shock  front  reflectivities  of  water 
and  nitromethane  could  be  very  important  from  both  fundamental  and  applied  view- 
points.  The  results  presented  in  this  report  suggest  the  possibility  that  explo¬ 
sive-related  chemistry  is  already  occurring  within  the  shock  front  thickness  of 
liquid  nitromethane,  and  that  the  effects  of  that  chemistry  can  be  observed  in 
the  optical  reflectivity  results.  Were  such  a  possibility  to  become  fact,  one 
would  have  a  new  tool  for  observing  and  controlling  detailed  explosive  properties 
in  terms  of  molecular  structure  and  composition.  From  a  fundamental  point  of 
view,  the  role  of  liquid  structure  and  composition  on  shock  front  structure  is 
fascinating  in  its  own  right. 


MULTIPLE  REFLECTION  THEORY 


In  our  previous  work  (ref  1)  the  effect  of  multiple  reflections  was 
neglected.  For  large  reflectivities  (e.g.,  20%  at  an  82°  optical  angle  of  inci¬ 
dence  in  water)  such  a  neglect  can  lead  to  difficulties.  Multiple  reflections 
are  considered  here  by  means  of  a  perturbative  approach. 


Consider  a  differential  region  of  varying  indexes  of  refraction  as  shown  in 
rigure  1* 

In  the  absence  of  multiple  reflections,  the  differential  reflected  amplitude 
is  given  by  ^ 


dEr  =  f  (Jr-)  (1  *  tan20) 


(1) 


where  0  is  the  angle  of  incidence  (measured  with  respect  to  the  normal  to  the 
surface  of  dn),  and  the  +,  -  signs  refer  to  perpendicular,  or  parallel  (to  the 
plane  of  incidence)  optical  polarizations,  respectively. 
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Figure  1.  Incoming  and  reflected  fields,  E  and  E  respectively,  impinging  on  a 
slab  of  differential  index  of  refraction  6n  ,  centered  at  Z,  , ,  and 
of  width  AZ  i+1  1+1 


Applying  equation  1  to  both  E  and  Er  in  figure  1  gives 

6n  • 

dEr(Zi+l}  =  2  [l  ±  tan20(Zi  +  +f-)  + 

i  6n 

-  2  (ITCZ^T]  I*  *  tan^e(Z1+1  +  f)]Er(Z1+1  +  f) 

Upon  expanding  terms  such  as  Er(Z  +  |£)  in  equation  2, 


(2) 


a  7  dE  d2E  9 

Vzi+1  +  =  VW  -arli+1(f)  +idF£l1+i(f)  + 


(3) 


gives  to  lowest  order 

dE  (E  -  E  ) 

dZ^ - 2n~ ~  ±  tan20)  (4) 
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Equation  4  simply  says  that,  to  lowest  order,  (E  -  E  )  replaces  E  in  equation  1 
when  rereflection  is  taken  into  account.  This  result  could,  of  course,  have  been 
guessed . 

Rewrite  equation  4  as 


[ffi  +  Is  £><1± 


(5) 


and  make  the  transformations 


Er  =  F(Z)ef(Z),  E  s  G(Z)ef(Z) 


(6) 


Substituting  equation  6  into  equation  5  gives 


dF  G  dn  . ,  ,  o„. 

dZ  =  2^  dZ  (1  1  tan  9) 


(7) 


if  f(Z)  is  taken  to  satisfy  the  differential  equation 


df 

dZ 


J-  (1  + 

2n  dZ  '  ~ 


tan20 ) 


(8) 


Thus,  when  (  f)  is  taken  to  be  the  reflected  amplitude  associated  with  unit 
amplitude  input  electric  field,  F  may  be  viewed  as  the  reflected  amplitude  asso¬ 
ciated  with  field  strength  G  in  a  space  where  rereflection  of  F  does  not  occur* 
That  is  the  beauty  of  the  transformation  given  by  equation  6;  one  can  solve  for  F 
exactly  as  in  the  simple  case  (ref  1)  where  multiple  reflections  are  not  impor¬ 
tant. 


For  a  shock  front  beginning  at  Z  =  0,  the  reflectivity,  R,  is  given  by 
2 


R  = 


E  (0) 

r 


E(0) 


_  if(°)  i2 
1 G(0) 1 


(9) 


Thus  F(0)  and  G(0)  determine  the  reflectivity. 

Poynting f s  theorem  for  the  reflection  problem  being  considered  here  can  be 
written  as 


CQ( |E(0) | 2  -  |Er(O)|2)sin0o  = 

=  C( | E | 2  -  |Er|2)sin9 


(10) 
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where  CQ  is  the  velocity  of  light  in  the  preshocked  medium. 


For 


v(Z)  =  C(Z)e2^Z) 


(ID 


Equation  10  becomes 


v  (|G(0)|2  -  | F(0) | 2)sin0  = 

o 

=  V( I G | 2  -  | F | 2)sin0 


(12) 


Applying  Snell’s  law,  nsinQ  =  n^inQ^  to  equation  12  gives 


F(Z), 

G(0)' 


2 


(13) 


One 

neglected 


now  proceeds  as  in  reference  1, 
by  defining  fractional  optical  phase 


where  multiple  reflections 

angles,  y  and  y  ,  by 
A  B 


were 


G(Z)  _  ,G(Z)_|  2lTlYA 
G(0)  ^(O)'6 


dF 

G(  0)  = 


(d| 


G(0) 


le2'1^ 


+  2iri| 


2iriy 


G(0) 


dy. 


Upon  rewriting  equation  13  to  lowest  order  in  I — —  I 

1 G(0) 1  * 


'G(O) 


=  { 


ne 


f(0) 


n  e 
o 


f(Z) 


■i  i’  -  m  i2i1/2u 


F  2  ^2 

+  i«f)i  i 2 


(14) 


(15) 


(16) 


and  substituting  that  along  with  equations  14  and  15  into  equation  7  gives 


=  rds.)  h  -  i£(o) 

l2nJ  11  lG(0) 


2^2 


(1  ±  tan20)cos{27r(y  -  y  )} 

A  B 


(17) 


The  solution  proceeds  by  solving  equations  8  and  17  along  with  an  equation  for 

kta  V* 
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Because  F  and  G  in  equation  7  appear  with  the  absence  rereflection,  writing 

down  equations  for  their  fractional  phase  angles  y  and  y  is  relatively  simple* ** 

A  B  r 


,  COS0 

d^B  =  X~  dz 


(18a) 


,  COS0 

dYA  =  —  dZ 


(18b) 


Integrating  equations  18  while  using  Snell's  law  and  X 
the  optical  wavelength  in  vacuum)  gives 


_1 

A  n  (with  \  being 
oo  oo 


<TA  -  V  -  - 


OO 


/q  {n2(Z')  -  n^sin2©^2  dZ' 


09) 


Again,  the  foot  of  the  shock  front  is  assumed  to  begin  at  Z  =  0,  and  again  our 
problem  is  solved  to  first  order  in  rereflection  by  simultaneously  solving  equa- 
tions  8,  17,  and  19. 


SOLUTIONS  WITH  REREFLECTION 


It  is  simple  to  show  that 


(1  ±  tan20)  =< 


n‘ 


n2  -  n2  sin20 
o  o 

n2  -  2n2  sin20 

_ o _ o 

n2  -  n2  sin20 
o  o 


so  that  integrating  equation  8  from  n^  to  n  gives 


n2  -  n2  sin20 

C+(Z)  -  --|  «n  {— - 2 - 


n2  -  n2  sin20 
o  o  o 


(20) 


(21) 


*  A  ^reflection  contribution  to  G(Z)  would,  in  principle,  complicate  the  phase 

ya  (Z)- 


5 


(22) 


f_(Z)  =  -  f+(Z)  +  tn(“) 


Equation  19  has  previously  (ref  2)  been  evaluated  for  the  special  case  of  a 
constant  gradient  in  the  index  of  refraction, 


n  =  nQ  +  (An)f  - 


(23) 


where  L  is  the  shock  front  thickness  and  (An)  is  the  total  change  in  index  of 

refraction  across  the  shock  front.  That  equation  19  result  for  the  special  case 
of  equation  23  is 


n2L 


(v  -  Y  )  - 2 - 

WA  V  y  (An) 


oo 


-  sin20  in 
o 


,  sin0  p  -  cos 0  + 

<"  +  1>2  V*-(-5r a) 


f 

/  sin0  p 

(m+1)  j 1  +  y 

k  -  {  xl°) 

1  m+1  J 

1  +  cos0 


(24) 


where 


m 


(25) 


has  been  used.  ra  is  the  most  convenient  parameter  to  use  for  the  numerical  inte¬ 
gration  of  the  right  hand  side  of  equation  17. 

As  Z  goes  from  Z  =  0  to  the  top  of  the  shock,  | -■  ■■  |  goes  from  to  0 

(reflection  ceases  at  the  top  of  the  shock).  Thus 


F(°)  -  j - ■)  =  An  {R1/2  +  /m  } 

1  +  \m  I 


(26) 


where  equation  9  has  been  used,  and  equation  17  becomes 


An  {R1/2  +  /l  +  R 

/I  -  R 


1 

2n 

o 


n  f(0) 

fn  dn  {  ~fC£j  1  n  ±  tan20)cos{  2tt(ya  -  Yg)} 
o  e 


(27) 
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Equation  27  should  be  compared  with  the  previous  (ref  1,  2)  result  where 

rereflection  and  F(Z)  on  the  right  hand  side  of  equation  12  are  neglected.  That 
previous  result  is 


R  =  sin  ["2“ 
o 


nf 

/  dn 
J  n 

o 


(1  ±  tan2  0  )  cos  {  2tt  (y  -  y  )}] 

A  Jd 


(28) 


In  the  limit  of  n  +  n^  and  R  -*■  0  equations  27  and  28  give  indentical  results. 

Let  Lj.  denote  the  right  hand  side  of  equation  27.  To  order  R3/2,  equation 
27  may  then  be  written  as 


R 


1/2 


(1 


rn 

+  t)  =  - 


which  becomes  (to  order  R  ) 


R2  +  ^R  -  j  =  0 


Solving  equation  30  for  R  to  order  I  4  then  gives 


R  =  I?  - 


±  3  •L± 


(29) 


(30) 


(31) 


Equation  31,  along  with  equations  21,  22,  and  24,  represents  the  final 
result  of  this  section.  To  the  algebraic  order  considered,  the  result  is  consis¬ 
tent  with  a  first  order  correction  for  multiple  reflections  (i.e.,  for  the  pre¬ 
sence  of  rereflection)  and  finite  reflectivity  (e.g.,  F(Z)  in  the  right  hand  side 
of  equation  12  taken  to  be  non-zero).  The  algebraic  accuracy  of  equation  31  is 
such  that  inclusion  of  the  next  higher  order  term  for  R  in  equation  29  is  equiva¬ 
lent  to  a  1%  reflectivity  calculational  error  when  R  =  0.2  (i.e.  20%  reflec¬ 
tion). 


It  is  a  fairly  simple  task  to  evaluate  I_j_  numerically  by  means  of  parti¬ 
tions  and  a  cyclic  computer  program.  Indeed,  a  19-partition  integration  program 
was  run  on  a  TRS-80  handheld  computer  with  each  calculation  point  taking  approxi¬ 
mately  5  minutes.  The  results  of  those  calculations  are  shown  and  discussed  in 
the  next  section.  The  calculational  error  associated  with  19  as  compared  39  pro¬ 
gram  steps  throughout  the  entire  range  of  experimental  interest  is  less  than  1%. 
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EXPERIMENTAL  AND  THEORETICAL  RESULTS 


The  experimental  data  points  shown  in  figures  2,  3,  and  4  are  the  results  of 
TaV?o^imental  program  at  the  Laboratory  for  Energetics  and  Detonation  (C.N.R.S. 
, a  3'  of  the  University  of  Poitiers.  That  experimental  program  is  supported 
by  the  Centre  National  de  la  Recherche  Scientifique  (C.N.R.S.)  of  France.  The 
theory  presented  in  the  previous  section  is  supported  by  the  U.S.  Army  through 
he  Li-House  Laboratory  Independent  Research  (ILIR)  Program  within  the  Large 
Caliber  Weapons  Systems  Laboratory  (LCL)  of  ARRADCOM. 


The  experimental  data  points  were  obtained  with  the  small  light-gas  canon 

^  3  t*16  University  °L  Poitiers.  A  stabilized  argon  laser  operating  at 

4145A  was  used  to  illuminate  the  shock  front,  and  an  8-mm  diameter,  10-namosecond 
rise  time,  silicon  photodiode  was  used  to  measure  the  calibrated  reflectivity. 

e  pressure  determination  (5.8  kbar  in  water  and  6.0  kbar  in  nitromethane)  was 
measured  to  within  5%,  while  the  reflectivity  record  was  read  to  within  10%. 

The  authors  have  previously  compared  (ref  1)  experiment  with  theory  (without 
multiple  reflections)  for  water.  With  multiple  reflections  (fig.  2) ,  8  <0.1 

fits  the  data  fairly  well,  while  a  somewhat  larger  value  of  g  appears1  to  be 

necessary  for  the  parallel  to  the  plane-of-incidence  optical  polarization.  In 
the  absence  of  multiple  reflections,  the  theory  fits  both  optical  polarizations 
equally  well  (g  <  0.1).  While  it  cannot  be  said  at  this  time  whether  the  differ¬ 
ence  in  agreement  between  theory  and  experiment  for  the  two  optical  polarizations 

in  water  has  a  real  physical  basis,  in  principle,  the  two  optical  polarizations 

should  not  quite  “see"  the  same  thing* 


Rahman  and.  Stillinger  (ref  3)  have  demonstrated  (by  computer  molecular 
dynamics)  a  finite  lifetime  for  shear  waves  in  liquid  water  (e.g.,  ~0.4  x  10“^ 2 
sec  for  a  shear  wavelength  of  ~13A).  The  implication  of  that  finite  lifetime  is 
the  initial  one-dimensional  strain's  becoming  hydrostatic  only  after  the  shock 
front  (i.e.,  after  the  shock  plateau)  has  been  reached.  The  0.4  x  10-12  sec 
lifetime  is  consistent  with  a  shock  front  rise  time  in  water  of  ~10-^2  sec  as  has 
independently  been  determined  (ref  4)  from  shock  polarization  theory  and  experi¬ 
ment.  Thus,  the  perpendicular  optical  polarization  case  (with  E  perpendicular  to 
the  direction  of  shock  propagation),  and  the  parallel  polarization  case  (where 

^S,,  component  of  E  parallel  to  the  direction  of  shock  propagation)  would 
not  see  identical  physical  phenomena. 


The  data  and  theory  (with  rereflection  included)  for  liquid  nitromethane  are 
s  own  in  figure  3.  The  systematic  loss  of  agreement  between  theory  and  experi- 
ment,  for  both  polarizations,  with  increasing  angle  of  incidence  (decreasing 
grazing  angle)  is  clear  and  graphic.  The  approximately  50%  difference  between 
theory  and  experiment  at  8  degrees  grazing  angle  is  well  beyond  theoretical  and 
experimental  uncertainty  and  error. 


The  optical  path  length  through  the  shock  front  is  given  by  (2L/cos0)  where 
L  is  the  shock  front  thickness.  The  systematic  worsening  of  agreement  between 
theory  and  experiment  suggests  the  possible  existence  of  an  optical  path  length 
effect  in  the  reflectivity  measurements.  In  turn,  an  optical  path  length  effect 
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REFLECTIVITY  (R) 


WATER,  5.8  KBAR 


Figure  2.  Experimental  and  theoretical  (from  equation  31)  reflectivities  for  mp 
0.0394.  8  =  n0L/A00.  Solid  lines  are  theory,  and  X  and  •  are  experi¬ 
mental  data  points.  The  subscripts  on  8  inidcate  the  optical  polari¬ 
zation.  For  liquid  water.  Grazing  angle  is  (0  -  90°) 
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REFLECTIVITY  (R) 


10-1 


IQ'2 


Figure  3. 


NITROMETHANE,  6.0  KBAR 

X  (PERPENDICULAR  OPTICAL  POLARIZATION) 
•  (PARALLEL  OPTICAL  POLARIZATION) 


5  10  15  20  25 

GRAZING  ANGLE  (DEGREES) 


0.2 


0.2 


Experimental  and  theoretical  (from  equation  31)  reflectivities  for  m^  = 
0.0394.  6  =  n0L/A0o*  Solid  lines  are  theory,  and  X  and  •  are  experi¬ 

mental  data  points.  The  subscripts  on  g  indicate  the  optical  polari¬ 
zation.  For  liquid  nitromethane .  Grazing  angle  is  ( 0O  -  90°) 
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REFLECTIVITY  (R) 


10- 


10-2 


Figure  4. 


X  and  •  are  experimental' data  points.  Solid  curves  are  the  effect  of 
"Thermal  Fluctuations"  added  to  the  g  =  0.1  curves  of  figure  3.  For 
liquid  nitromethane.  Grazing  angle  is  (0  -  90°) 
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suggests  the  possible  existence  of  an  optical-scattering-related  fluctuation 
phenomenon. 

Since  nitrome thane  is  an  explosive,  it  is  conceivable  that  the  reflectivity 
data  are  demonstrating  the  existence  of  chemically  related  thermal  fluctuations 
within  the  shock  front.  Toops  has  reported  (ref  5)  on  the  temperature  dependence 
of  the  index  of  refraction  for  liquid  nitromethane .  His  observations  give 

dn  ,  -4  n  -1 

dT  x  10  c  (32) 


so  that  a  (AT)  >  0  fluctuation  will  decrease 
contribution  at  the  fluctuation  site. 


the  measured  index  of  refraction 


Consider  a  reflectivity,  R,  approximately  given  by 


R1/2 


(An) 

—  (l±tan20) 


and  a  fluctuation  effect  on  (An)^  given  by 


(33) 


(An)f  =  (An)£o)  (1  +  ct(T)h)  (34) 

where  h  =  2Lcos  *0  is  the  optical  path  length.  If  (AR)  is  taken  as  the  change 
( ot h / ^ ^ C ^^hen" Y  C3USed  by  the  thermal  fluctuations  characterized  by  a,  and  if 


(AR) 


9l 


(AR)02 


1  +  tan^Qj 
^l"  +  tanz02^ 


COS02 

COS01 


(35) 


If  (AR)e  is  taken  as  the  difference  between  experimental  and  theoretical 
values  for  B  u  =  0.1  at  0  =  82°,  then 


i(AR)ooO  =  -  5.5  x  10 


-2 


1 1 ( AR) 82O  =  “  4.4  x  10 


-2 


(36a) 

(36b) 


i^fiture11?  eq^ftl0"S  36a  aad36b  into  equation  35  yields  the  solid  curves  shown 
g  re  4.  Clearly  a  path  length  dependent  subtractive  contribution  to  An  by 
means  of  equations  33  through  35  serves  to  facilitate  the  agreement  between  the¬ 
ory  and  experiment. 
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DISCUSSION  AND  CONCLUSIONS 


The  most  striking  aspect  of  the  theoretical  and  experimental  results  pre¬ 
sented  in  this  report  is  the  distinct  difference  between  theoretical  and  experi¬ 
mental  agreement  for  liquid  water  and  liquid  nitromethane .  While  the  agreement 
for  water  is  not  perfect,  the  agreement  for  nitromethane  is  clearly  lacking  in 
the  absence  of  some  angle  of  incidence  dependent  effect  such  as  that  proposed  by 
equations  33  and  34.  Liquid  nitromethane  is  an  explosive  and  it  seems  reasonable 
to  postulate  explosive-related  chemistry  (although  6  kbar  is  approximately  an 
order  of  magnitude  less  than  room-temperature  shock  initiation  pressure)  as  the 
mechanism  for  the  a(T)h  term  of  equation  34. 


If  the  experimental-theoretical  discrepancy  shown  in  figure  3  can  be  attri¬ 
buted  to  explosive-related  chemical  effects,  as  seems  very  possible,  the  authors 
believe  that  the  attribution  would  be  a  first  for  the  physics  of  explosive  (ener¬ 
getic)  media.  If  explosive— related  chemistry  is  the  mechanism  for  proposed  ther— 
mal  fluctuations,  then  the  experimental-theoretical  discrepancy  should  increase 
with  increasing  shock  pressure.  Time  and  further  experiments  will  tell. 

mf  =  0.0394  was  used  for  both  water  at  5.8  kbar,  and  nitromethane  at  6.0 
kbar.  For  liquid  nitromethane  at  20°C,  from  Toops  (ref  5)  and  Hardesty  (ref  6), 

An  =  0.338  (p  -  p^)  (37) 

where  p  is  mass  density.  From  Enig  and  Petrone  (ref  7)  p  (6  kbar)  =  1.304  g  cm"3 
so  that  (An)f  =  0.0544.  The  preshock  value  of  nQ  =  1.385  then  gives  mf  =  0.0393. 

For  water,  from  Zeldovich,  et  al  (ref  8), 

An  =  0.334  (p  -  1) 

“3 

From  Ap  =  0.1575  g  cm  (an  extrapolation  of  the  Rice  and  Walsh  (ref  9)  data)  and 
n?  -  i*334  one  finds  m^  =  0.0394.  As  a  0.3%  change  in  m^  can  barely  be  seen  on  a 
plot  of  theoretical  results,  the  same  m^  value  was  used  in  characterizing  the 
index  of  refraction  change  in  both  liquids .  Similarly  a  10°C  shock-induced  tem¬ 
perature  change  has  a  negligible  effect  on  the  theoretical  predictions  of  equa¬ 
tion  32. 
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